Nonlinear harmonic generation and devices in doubly-resonant Kerr cavities 
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We describea theoretical analysis of the nonlinear dynamics of third-harmonic generation (cj — > 
3oo) via Kerr (x^) nonlinearities in a resonant cavity with resonances at both oo and 3a;. Such a 
doubly resonant cavity greatly reduces the required power for efficient harmonic generation, by a 
factor of ~ V/Q 2 where V is the modal volume and Q is the lifetime, and can even exhibit 100% 
harmonic conversion efficiency at a critical input power. However, we show that it also exhibits a rich 
variety of nonlinear dynamics, such as mult ist able solutions and long-period limit cycles. We describe 
how to compensate for self/cross-phase modulation (which otherwise shifts the cavity frequencies 
out of resonance), and how to excite the different stable solutions (and especially the high-efficiency 
solutions) by specially modulated input pulses. 

PACS numbers: 42.65.Ky, 42. 60. Da, 42.65.Sf, 42.65.Jx 



I. INTRODUCTION 



In this paper, we describe how 100% third-harmonic 
conversion can occur in doubly-resonant optical cavities 
with Kerr nonlinearities, even when dynamical stabil- 
ity and self-phase modulation (which can drive the cavi- 
ties out of resonance) are included (extending our earlier 
work [1]), and describe the initial conditions required to 
excite these efficient solutions. In particular, we show 
that such doubly-resonant nonlinear optical systems can 
display a rich variety of dynamical behaviors, including 
multistability (different steady states excited by vary- 
ing initial conditions, a richer version of the bistable 
phenomenon observed in single- mode cavities [2]), long- 
period limit cycles (similar to the "self-pulsing" observed 
for second-harmonic generation [3(]), and transitions in 
the stability and multiplicity of solutions as the parame- 
ters are varied. One reason for studying such doubly res- 
onant cavities was the fact that they lower the power re- 
quirements for nonlinear devices [1], and in particular for 
third harmonic conversion, compared to singly-resonant 
cavities or nonresonant structures MU5|, [y, |7|, |8|, |9|, LLO, LUU, 
HI HE HI HE HE Hi HI, HI, HE EOT An appreciation 
and understanding of these behaviors is important to de- 
sign efficient harmonic converters (the main focus of this 
paper), but it also opens the possibility of new types of 
devices enabled by other aspects of the nonlinear dynam- 
ics. 

In a Kerr (x^) medium, there is a change in the refrac- 
tive index proportional to the square of the electric field; 
for an oscillating field at a frequency a;, this results in a 
shift in the index at the same frequency (self-phase modu- 
lation, SPM), generation of power at the third-harmonic 
frequency 3a;, and also other effects when multiple fre- 
quencies are present [cross-phase modulation (XPM) and 
four- wave mixing (FWM)] [22j]. When the field is con- 
fined in a cavity, restricting to a small modal volume V 
for a long time given by the quality factor Q (a lifetime 
in units of the optical period) [23], such nonlinear effects 
are enhanced by both the increased field strength for the 



same input power and by the frequency sensitivity inher- 
ent in resonant effects (since the fractional bandwidth is 
1/Q). This enhancement is exploited, for example, in 
nonlinear harmonic and sum-frequency generation, most 
commonly for x^ effects where the change in index is 
proportional to the electric field (which requires a non- 
centrosymmetric material) [22]. One can further enhance 
harmonic generation by using a cavity with two resonant 
modes, one at the source frequency and one at the har- 
monic frequency [1 HI HE EE HE HE HE HO, H2 . In this 
case, one must also take into account a nonlinear down- 
conversion process that competes with harmonic genera- 
tion (EHEUli, but it turns out to be theoretically pos- 
sible to obtain 100% harmonic conversion for either x^ 
(a; — > 2u) or x^ —> 3a; ) nonlinearities at a specific 
"critical" input power P cr it (both in a one-dimensional 
model of propagating waves for x^ nonlinearities [32[ 
and also in a more general coupled-mode model for either 
^( 2 ) or x^ nonlinearities [l|). In particular, we studied 
the harmonic-generation and downconversion processes 
in a broad class of model systems depicted in Fig. [TJ 
a single input channel (e.g. a waveguide) is coupled to 
a nonlinear cavity with two resonant frequencies, where 
both reflected and harmonic fields are emitted back into 
the input channel. In this case, we predicted 100% har- 
monic generation at a critical power P cr it proportional 
to V/Q 3 for x (2) and V/Q 2 for x (3) Q- However, we 
only looked at the steady-state solution of the system and 
not its dynamics or stability. Moreover, in the x^ case 
there can also be an SPM/XPM effect that shifts the cav- 
ity frequencies out of resonance and spoils the harmonic- 
generation effect. In this paper, we consider both of these 
effects, describe how to compensate for SPM/XPM, and 
demonstrate the different regimes of stability in such x^ 
doubly resonant systems. We show that the parameters 
and the initial conditions must be chosen within certain 
regimes to obtain a stable steady state with high conver- 
sion efficiency. 

In other regimes, we demonstrate radically different 
behaviors: not only low-efficiency steady states, but also 
limit-cycle solutions where the efficiency oscillates slowly 
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FIG. 1: Top: Schematic of general scheme for third-harmonic 
generation, and dynamical variables for coupled-mode equa- 
tions: a single input/output channel (with incoming/outgoing 
field amplitudes s±) is coupled to a resonant cavity with two 
modes at frequencies u)i and 3ui (and corresponding ampli- 
tudes a\ and 0-3). The two resonant modes are nonlinear ly 
coupled by a Kerr (x^) nonlinearity. Bottom: An exam- 
ple realization in one dimension, using a semi-infinite 
quarter- wave stack of dielectric layers with a doubled-layer de- 
fect (resonant cavity) that is coupled to incident plane waves; 
the electric field of a steady-state 3cji solution is shown as 
blue/white/red for negative/zero/positive. 



with a repetition period of many thousands of optical cy- 
cles. With infrared light, these limit cycles form a kind 
of optical oscillator/clock with a period in the hundreds 
of GHz or THz (and possibly lower, depending on the 
cavity parameters). Previously, limit-cycle/self-pulsing 
behaviors have been observed in a number of other non- 
linear optical systems, such as: doubly-resonant x^ cav- 
ities coupled by second-harmonic generation [3[ ; bistable 
multimode Kerr cavities with time-delayed nonlineari- 
ties (33t ; nonresonant distributed feedback in Bragg grat- 
ings [3J]; and a number of nonlinear lasing devices [351 ]. 
However, the system considered in this work seems un- 
usually simple, especially among x^ systems, in that 
it only requires two modes and an instantaneous Kerr 
nonlinearity, with a constant-frequency input source, to 
attain self-pulsing, and partly as a consequence of this 
simplicity the precise self-pulsing solution is quite insen- 
sitive to the initial conditions. In other nonlinear optical 
systems where self-pulsing was observed, other authors 
have also observed chaotic solutions in certain regimes. 
Here, we did not observe chaos for any of the parame- 
ter regimes considered, where the input was a constant- 
frequency source, but it is possible that chaotic solutions 
may be excited by an appropriate pulsed input as in the 
X (2) case 

Another interesting phenomenon that can occur in 
nonlinear systems is mult ist ability, where there are mul- 
tiple possible steady-state solutions that one can switch 
among by varying the initial conditions. In Kerr (x^) 
media, an important example of this phenomenon is 
bistable transmission through nonlinear cavities: for 
transmission through a single-mode cavity, output can 
switch discontinuously between a high-transmission and 
a low-transmission state in a hysteresis effect that re- 
sults from SPM [2]. For example, if one turns on the 
power gradually from zero the system stays in the low- 



transmission state, but if the power is increased fur- 
ther and then decreased to the original level, the sys- 
tem can be switched to the high-transmission state. 
This effect, which has been observed experimentally (36[, 
can be used for all-optical logic, switching, rectification, 
and many other functions [2]. In a cavity with mul- 
tiple closely-spaced resonances, where the nonlinearity 
is strong enough to shift one cavity mode's frequency 
to another's, the same SPM phenomenon can lead to 
more than two stable solutions [37| • Here, we demon- 
strate a much richer variety of multistable phenomena in 
the doubly-resonant case for widely-separated cavity fre- 
quencies coupled by harmonic generation in addition to 
SPM — not only can there be more than two stable states, 
but the transitions between them can exhibit complicated 
oscillatory behaviors as the initial conditions are varied, 
and there are also Hopf bifurcations into self-pulsing so- 
lutions. 

The remaining part of the paper is structured as fol- 
lows. In Sec. HI1 we review the theoretical description of 
harmonic generation in a doubly resonant cavity coupled 
to input /output waveguides, based on temporal coupled- 
mode theory, and demonstrate the possibility of 100% 
harmonic conversion. We also discuss how to compen- 
sate for frequency shifting due to SPM and XPM by pre- 
shifting the cavity frequencies. In Sec. we analyze 
the stability of this 100%-emciency solution, and demon- 
strate the different regimes of stable operation that are 
achieved in practice starting from that theoretical ini- 
tial condition. We also present bifurcation diagrams that 
show how the stable and unstable solutions evolve as the 
parameters vary. Finally, in Sec. IIV[ we consider how 
to excite these high-efficiency solutions in practice, by 
examining the effect of varying initial conditions and un- 
certainties in the cavity parameters. In particular, we 
demonstrate the multistable phenomena exhibited as the 
initial conditions are varied. We close with some con- 
cluding remarks, discussing the many potential directions 
for future work that are revealed by the phenomena de- 
scribed here. 



II. 100% HARMONIC CONVERSION IN 
DOUBLY-RESONANT CAVITIES 

In this section, we describe the basic theory of fre- 
quency conversion in doubly-resonant cavities with x^ 
nonlinear it ies, including the undesirable self- and cross- 
phase modulation effects, and explain the existence of a 
solution with 100% harmonic conversion (without con- 
sidering stability). 

Consider a waveguide coupled to a doubly resonant 
cavity with two resonant frequencies o;f av = uj\ and 
c^3 av = ujs = 3uji (below, we will shift cc^ av to differ 
slightly from c^), and corresponding lifetimes T\ and T3 
describing their radiation rates into the waveguide (or 
quality factors Qk — UkTk/2). In addition, these modes 
are coupled to one another via the Kerr nonlinearity. Be- 
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cause all of these couplings are weak, any such system 
(regardless of the specific geometry), can be accurately 
described by temporal coupled-mode theory, in which the 
system is modeled as a set of coupled ordinary differ- 
ential equations representing the amplitudes of the dif- 
ferent modes, with coupling constants and frequencies 
determined by the specific geometry p], [38[ . In particu- 
lar, the coupled- mode equations for this particular class 
of geometries were derived in Ref. 1 along with explicit 
equations for the coupling coefficients in a particular ge- 
ometry. The degrees of freedom are the field amplitude 
dk of the kth cavity mode (normalized so that \dk\ 2 is the 
corresponding energy) and the field amplitude Sk± of the 
incoming (+) and outgoing (— ) waveguide modes at uou 
(normalized so that \su±\ 2 is the corresponding power), 
as depicted schematically in Fig. [TJ These field ampli- 
tudes are coupled by the following equations (assuming 
that there is only input at uj\, i.e. s 3 + = 0): 
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As explained in Ref. 1, the a and /3 coefficients are 
geometry /material-dependent constants that express the 
strength of various nonlinear effects for the given modes. 
The dij terms describe self- and cross-phase modulation 
effects: they clearly give rise to effective frequency shifts 
in the two modes. The term characterize the energy 
transfer between the modes: the (3 3 term describes fre- 
quency up-conversion and the _3\ term describes down- 
conversion. As shown in Ref. they are related to one 
another via conservation of energy uj\fi\ = uj 3 P^ and 
all of the nonlinear coefficients scale inversely with the 
modal volume V. 

There are three different ol^ parameters (two SPM co- 
efficients an and 0^33 and one XPM coefficient ol\ 3 — 
a 3 i). All three values are different, in general, but are 
determined by similar integrals of the field patterns, pro- 
duce similar frequency-shifting phenomena, and all scale 
as 1/V. Therefore, in order to limit the parameter space 
analyzed in this paper, we consider the simplified case 
where all three frequency-shifting terms have the same 
strength ol^ = a. 

One can also include various losses, e.g. linear losses 
correspond to a complex uj\ and/or co 3 , and nonlin- 
ear two-photon absorption corresponds to a complex a. 
As discussed in the concluding remarks, however, we 
have found that such considerations do not qualitatively 
change the results (only reducing the efficiency some- 
what, as long as the losses are not too big compared to 
the radiative lifetimes r), and so in this manuscript we 
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FIG. 2: Steady-state em(fi§J§$-| of third-harmonic generation 
(solid red line) from Ref. [0, for a = (no self-phase mod- 
ulation), as a function of input power |si+| 2 scaled by the 
Kerr coefficient 712 = /Ace. The reflected power at the 
incident frequency uji is shown as a dashed black line. There 
is a critical power where the efficiency of harmonic genera- 
tion is 100. The parameters used in this plot are Qi = 1000, 
Q 3 = 3000, Pi = (4.55985 - 0.7244i) x 10" 5 in dimensionless 
units of x (3) /^ e - 



restrict ourselves to the idealized lossless case. 

Figure [2] shows the steady-state conversion efficiency 
(| S3- | 2 /|si + | 2 ) versus input power of light that is inci- 
dent on the cavity at (jf™, for the same parameter regime 
in Ref. 1 (i.e. assuming negligible self- and cross-phase 
modulation so that a = 0), and not considering the sta- 
bility of the steady state. As shown by the solid red 
curve, as one increases the input power, the efficiency in- 
creases, peaking at 100% conversion for a critical input 
power P cr it = l^i"*!, where 



5 i+ l 



\ui(3l\ 2 T?T 3 



1/4 



(3) 



The efficiency decreases if the power is either too low 
(in the linear regime) or too high (dominated by down- 
conversion). P cr it scales as V/Q 2 , so one can in principle 
obtain very low-power efficient harmonic conversion by 
increasing Q and/or decreasing V [1]. Including absorp- 
tion or other losses decreases the peak efficiency, but does 
not otherwise qualitatively change this solution [l|. 

There are two effects that we did not previously ana- 
lyze in detail, however, which can degrade this promising 
solution: nonlinear frequency shifts and instability. Here, 
we first consider frequency shifts, which arise whenever 
« / 0, and consider stability in the next section. The 
problem with the a terms is that efficient harmonic con- 
version depends on the cavities being tuned to harmonic 
frequencies uo 3 = 3cji; a nonlinear shift in the cavity 
frequencies due to self/cross-phase modulation will spoil 
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FIG. 3: Shift in the resonant Ire^ilemsy 
input power, due to self /cross-phase modulation. (There is an 
identical shift in w 3 NL .) If the cavity is designed so that the 
linear (Pi n — > 0) frequencies are harmonics, the nonlinearity 
pushes the system out of resonance (lower blue line) as the 
power increases to the critical power for 100% efficiency. This 
is corrected by pre-shifting the cavity frequencies (upper green 
line) so that the nonlinear frequency shift pushes the modes 
into resonance at P C rit- 



this resonance. In principle, there is a straightforward so- 
lution to this problem, as depicted in Fig. 03 Originally 
(for a = 0), the cavity was designed to have the frequency 
uji in the linear regime, but with the effective cav- 

ity frequency c^ L (including self/cross-phase modulation 
terms) is shifted away from the design frequency as shown 
by the blue line. Instead, we can simply design the linear 
cavity to have a frequency c<;5; av slightly different from the 
operating frequency uji, so that self/cross-phase modula- 
tion shifts L exactly to u\ at the critical input power, 
as depicted by the green line in Fig. [3j Exactly the same 
strategy is used for by pre-shifting uo^ w . 
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Then, we substitute these critical fields into the coupled- 
mode equations for a 7^ 0, and solve for the new cavity 
frequencies o;£ av so as to cancel the a terms and make 
the a£ nt solutions still valid. This yields the following 



transformation of the cavity frequencies: 



More precisely, to compute the required amount of 
pre-shifting, we examine the coupled-mode equations 
Eqs. (QJEJ). First, we solve for the critical power P cr it 
assuming a = 0, as in Ref. 1, and obtain the correspond- 
ing critical cavity fields a£ rit : 



UJl 



1 - an|a; rit | 2 - a 13 \af^ 
lj 3 



1 - a 13 \af 



(6) 
(?) 



By inspection, when substituted into Eqs. (pj{2]) at the 
critical power, these yield the same steady-state solution 
as for a = 0. (There are two other appearances of uji 
and cos in the coupled-mode equations, in the /3 k terms, 
but we need not change these frequencies because that is 
a higher-order effect, and the derivation of the coupled- 
mode equations considered only first-order terms in x^ •) 
The nonlinear dynamics turn out to depend only on 
four dimensionless parameters: ts/ti = Q3/3Q1, an/ft, 
0^33/^1, and ais//3i = a 3 i / (3\. The overall scale of Q, 
a, etcetera, merely determines the absolute scale for the 
power requirements: it is clear from the equations that 
multiplying all a and (3 coefficients by an overall con- 
stant K can be compensated by dividing all a and s 
amplitudes by \[K [which happens automatically for s 
at the critical power by Eq. (|3J)]; the case of scaling 
Ti ,3 by an overall constant is more subtle and is consid- 
ered below. As mentioned above, for simplicity we take 
an = <^33 = <^13 = <^3i = ol. Therefore, in the subse- 
quent sections we will analyze the asymptotic efficiency 
as a function of T3 / n and a /Pi. 

So far, we have found a steady-state solution to the 
coupled- mode equations, including self/cross-phase mod- 
ulation, that achieves 100% third-harmonic conversion. 
In the following sections, we consider under what condi- 
tions this solution is stable, what other stable solutions 
exist, and for what initial conditions the high-efficiency 
solution is excited. 

To understand the dynamics and stability of the non- 
linear coupled- mode equations, we apply the standard 
technique of identifying fixed points of the equations and 
analyzing the stability of the linearized equations around 
each fixed point [301 ] . 

By a "fixed point," we mean a steady-state solution 



corresponding to an input frequency uji (si + 



-%UJ\t 



and hence ai(t) = Aie~ iuJlt and a 3 (t) = A 3 e- i3uJlt for 
some unknown constants A\ and A%. [An input fre- 
quency uji can also generate higher harmonics, such as 
9cji or 5cji, but these are negligible: both because they 
are higher-order effects (~ [x^] 2 ? an d all such terms 
were dropped in deriving the coupled- mode equations), 
and because we assume there is no resonant mode present 
at those frequencies.] By substituting this steady-state 
form into Eqs. (p]-[2]), one obtains two coupled polynomial 
equations whose roots are the fixed points. We already 
know one of the fixed points from the previous section, 
the 100% efficiency solution, but to fully characterize the 
system one would like to know all of the fixed points 
(both stable and unstable). We solved these polynomial 
equations using Mathematica, which is able to compute 
all of the roots, but some transformations were required 
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to put the equations into a solvable form. In particu- 
lar, we eliminated the complex conjugations by writing 
Ak = r^e % ^ k and assuming (without loss of generality) 
that si+ is real. Multiplying Eq. (PQ) by e~ i(t>1 and Eq. p|) 
by e - * 3 ^ 1 allows us to simply solve the system for e - *^ 1 
and e*(03-30i)^ Requiring the magnitude of these two 
quantities to be unity yields two polynomials in x = r\ 
and y = r|, which Mathematica can handle. The result- 
ing polynomial is of an artificially high degree, resulting 
in spurious roots, but the physical solutions are easily 
identified by the fact that x and y must be real and non- 
negative. (We should also note that this root-finding 
process is highly sensitive to roundoff error [40[, inde- 
pendent of the physical stability of the solutions, but we 
dealt with that problem by employing 50 decimal places 
of precision.) 

As mentioned above, the dynamics are independent of 
the overall scale of 71,3, and depend only on the ratio 
T3/71. This can be seen from the equations for Ai^, 
in which the uj\^ oscillation has been removed. In these 
equations, if we multiply t\ and T3 by an overall constant 
factor K, after some algebra it can be shown that the 
A\$ equations are invariant if we rescale A\ — ► Ai/y/K, 
As — > As/\ /r K^ rescale time t — > Kt, and rescale the 
input si_|_ — » s\+/K [which happens automatically for 
the critical power by Eq. (j3])]. Note also that the con- 
version efficiency |<S3_/si + | 2 = (2/ts) \As/si+ | 2 is also 
invariant under this rescaling by K. That is, the powers 
and the timescales of the dynamics change if you change 
the lifetimes, unsurprisingly, but the steady states, sta- 
bility, etcetera (as investigated in the next section) are 
unaltered. 
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FIG. 4: Phase diagram of the nonlinear dynamics of the 
doubly-resonant nonlinear harmonic generation system from 
Fig. [T]as a function of the relative cavity lifetimes {t^/ti — 
3Q3/Q1) and the relative strength of SPM/XPM vs. har- 
monic generation (a//3±) for input power equal to the critical 
power for 100% efficiency. For T3 < n there is always one 
stable 100%-efficiency solution, and for nonzero a the system 
may have additional stable solutions. For T3 > n the 100%- 
efficiency solution becomes unstable, but there are limit cycles 
and lower-efficiency stable solutions. Various typical points 
A-G in each region are labeled for reference in the subsequent 
figures. 



III. STABILITY AND DYNAMICS 

Given the steady-state solutions (the roots), their sta- 
bility is determined by linearizing the original equations 
around these points to a first-order linear equation of 
the form dx/dt = Ax; a stable solution is one for which 
the eigenvalues of A have negative real parts (leading 
to solutions that decay exponentially towards the fixed 
point) [H. 

The results of this fixed-point and stability analysis are 
shown in Fig. Has a "phase diagram" of the system as 
a function of the relative lifetimes T3/71 = 3Q3/Q1 and 
the relative strength of self-phase- modulation vs. four- 
wave mixing a/f3\. Our original 100%-efficiency solu- 
tion is always present, but is only stable for T3 < t\ 
and becomes unstable for T3 > r±. The transition point, 
73 = Ti, corresponds to equal energy |ai| 2 = \as\ 2 in 
the fundamental and harmonic mode at the critical input 
power. The unstable region corresponds to 1 0,3 1 2 > |ai| 2 
(and the down-conversion term is stronger than the up- 
conversion term) — intuitively, this solution is unstable 
because, if any perturbation causes the energy in the 
harmonic mode to decrease, there is not enough pump- 
ing from up-conversion to bring it back to the 100%- 



efficiency solution. Conversely, in the stable |a-3 1 2 < |ai| 2 
(73 < T\) regime, the higher-energy fundamental mode 
is being directly pumped by the input and can recover 
from perturbations. Furthermore, as a//3i increases, ad- 
ditional lower-efficiency stable solutions are introduced, 
resulting in regimes with two (doubly stable) and three 
(triply stable) stable fixed points. These different regimes 
are explored in more detail via bifurcation diagrams be- 
low, and the excitation of the different stable solutions is 
considered in the next section. 

For T3 > ti, the 100%-efficiency solution is unstable, 
but there are lower-efficiency steady-state solutions and 
also another interesting phenomenon: limit cycles. A 
limit cycle is a stable oscillating-efficiency solution, one 
example of which (corresponding to point D in Fig. 
is plotted as a function time in Fig. [5l (In general, the 
existence of limit cycles is difficult to establish analyti- 
cally (39[, but the phenomenon is clear in the numerical 
solutions as a periodic oscillation insensitive to the ini- 
tial conditions). In fact, as we shall see below, these 
limit cycles result from a "Hopf bifurcation," which is a 
transition from a stable fixed point to an unstable fixed 
point and a limit cycle [4l|. In this example at point D, 
the efficiency oscillates between roughly 66% and nearly 
100%, with a period of several thousand optical cycles. 
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FIG. 5: An example of a limit-cycle solution, with a period- 
ically oscillating harmonic-generation efficiency as a function 
of time, corresponding to point D in Fig. 2] Perturbations in 
the initial conditions produce only phase shifts in the asymp- 
totic cycle. Here, the limit cycle has a period of around 3 x 10 4 
optical cycles. Inset: Square of Fourier amplitudes (arbitary 
units) for each harmonic component of the limit cycle in the 
Fourier- series expansion of the \As\. 




FIG. 6: Bifurcation diagram showing the harmonic- 
generation efficiency of the stable (solid red lines) and un- 
stable (dashed blue lines) steady-state solutions as a function 
oia/01 for a fixed t$/t\ =0.7, corresponding to the line ACF 
in Fig. [4] (see inset). The input power is the critical power 
Pcrit, so there is always a 100%-efficiency stable solution, but 
as a/ Pi increases new stable and unstable solutions appear 
at lower efficiencies. 



As a consequence of the time sacling described in the 
last paragraph of the previous section, the period of such 
limit cycles is proportional to the r's. If the frequency 
(j i were 1.55 /xm, for a Qi of 500 optical cycles, this limit 
cycle would have a frequency of around 70 GHz, forming 
an interesting type of optical "clock" or oscillator. Fur- 
thermore, the oscillation is not sinusoidal and contains 
several higher harmonics as shown in the inset of Fig. [5j 
the dominant frequency component in this case is the 
fourth harmonic (~ 280 GHz), but different points in the 
phase diagram yield limit cycles with different balances 
of Fourier components. 

To better understand the phase diagram of Fig. 2J it is 
useful to plot the efficiencies of both the stable and unsta- 
ble solutions as a function of various parameters. Several 
of these bifurcation diagrams (in which new fixed points 
typically appear in stable/unstable pairs) are shown in 
Figs. EHH1 To begin with, Figs. [6] and [71 correspond to 
lines connecting the labeled points ACF, BCD, and ECG, 
respectively, in Fig. [H showing how the stability changes 
as a function of a/ Pi and t^/t\. Figure [6] shows how first 
one then two new stable fixed points appear as a/ Pi is 
increased, one approaching zero efficiency and the other 
closer to 50%. Along with these two stable solutions ap- 
pear two unstable solutions (dashed lines). (A similar 
looking plot, albeit inverted, can be found in Ref. [3?] for 
SPM-coupled closely-spaced resonances.) In particular, 
the fact that one of the unstable solutions approaches 
the 100%-efficiency stable solution causes the latter to 



have a smaller and smaller basin of attraction as a/ pi 
increases, making it harder to excite as described in the 
next section. The next two plots, in Fig. [7j both show 
the solutions with respect to changes in T3 / t\ at two dif- 
ferent values of a/ pi . They demonstrate that at t\ = T3 , 
a Hopf bifurcation occurs where the 100%-efficiency solu- 
tion becomes unstable for T3 > r± and limit cycles appear, 
intuitively seeming to "bounce between" the two nearby 
unstable fixed points. (The actual phase space is higher 
dimensional, however, so the limit cycles are not con- 
strained to lie strictly between the efficiencies of the two 
unstable solutions.) It is worth to note that the remain- 
ing nonzero-efficiency stable solution (which appears at 
a nonzero T3/71) becomes less efficient as T3/71 increases. 

The above analysis and results were for the steady- 
state-solutions when operating at the critical input power 
to obtain a 100%-efficiency solution. However, one can, 
of course, operate with a different input power — although 
no other input power will yield a 100%-efficient steady- 
state solution, different input powers may still be useful 
because (as noted above and in the next section) the 
100%-efficiency solution may be unstable or practically 
unattainable. Figure [SJleft) is the bifurcation diagram 
with respect to the input power Pi n /P CT it at fixed a/ pi 
and fixed T3/71, corresponding to point C in Fig.JH This 
power bifurcation diagram displays a number of interest- 
ing features, with the steady-state solutions transitioning 
several times from stable to unstable and vice versa. As 
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FIG. 7: Bifurcation diagram showing the harmonic- generation efficiency of the stable (solid red lines) and unstable (dashed 
blue lines) steady-state solutions as a function of r^/n for a fixed a/ Pi = 3 (left) or = 8 (right), corresponding to the lines 
BCD or EFG, respectively, in Fig. |4] (see insets). The input power is the critical power P cr it, so there is always a 100%-efficiency 
steady-state solution, but it becomes unstable for T3 > n (a Hopf bifurcation leading to limit cycles as in Fig. [5}. 



we will see in the next section, the stability transitions 
in the uppermost branch are actually supercritical (re- 
versible) Hopf bifurcations to/from limit cycles. Near the 
critical power, there is only a small region of stability of 
the near- 100%-efficiency solution, as shown in the inset 
of Fig. [SK left). In contrast, the lower-efficiency stable so- 
lutions have much larger stable regions of the curve while 
still maintaining efficiencies greater than 70% at low pow- 
ers comparable to P cr i t ~ V/Q 2 , which suggests that they 
may be attractive regimes for practical operation when 
a /f3i is not small. This is further explored in the next 
section, and also by Fig. [Htright) which shows the bifur- 
cation diagram along the line ACF in Fig. 0] [similar to 
Fig. [6], but at 135% of the critical input power. For this 
higher power, the system becomes at most doubly sta- 
ble as a/ /3i is increased, and the higher-efficiency stable 
solution becomes surprisingly close to 100% as a /Pi — > 0. 



IV. EXCITING HIGH-EFFICIENCY 
SOLUTIONS 

One remaining concern in any multistable system is 
how to excite the desired solution — depending on the ini- 
tial conditions, the system may fall into different stable 
solutions, and simply turning on the source at the criti- 
cal input power may result in an undesired low-efficiency 
solution. If a/P is small enough, of course, then from 
Fig. [H the high-efficiency solution is the only stable solu- 
tion and the system must inevitably end up in this state 
no matter how the critical power is turned on. Many in- 
teresting physical systems will correspond to larger val- 
ues of a/P, however [l|, and in this case the excitation 
problem is complicated by the existence of other stable 
solutions. Moreover, the basins of attraction of each sta- 
ble solution may be very complicated in the phase space, 



as illustrated by Fig. [9l where varying the initial cavity 
amplitudes Ai^ from the 100%-efficiency solution causes 
the steady state to oscillate in a complicated way between 
the three stable solutions (at point C in Fig. 2]). We have 
investigated several solutions to this excitation problem, 
and found an "adiabatic" excitation technique that reli- 
ably produces the high-efficiency solution without unrea- 
sonable sensitivity to the precise excitation conditions. 

First, we considered a simple technique similar to the 
one described in Ref. [2| for exciting different solutions of 
a bistable filter: as shown in Fig.[10l we "turn on" the in- 
put power by superimposing a gradual exponential turn- 
on (asymptoting to Pi = P cr it) with a Gaussian pulse of 
amplitude Po and width ST. The function of the initial 
pulse is to "kick" the system into the desired stable solu- 
tion. We computed the eventual steady-state efficiency 
(after all transient effects have disappeared) as a function 
of the pulse amplitude Po at point C in Fig. 0J where 
there are three stable solutions. The results are shown in 
Fig. [TIT left), and indeed we see that all three stable solu- 
tions from point C in Fig. [6j one at near-zero efficiency, 
one at around 47% efficiency, and one at 100% efficiency. 
Unfortunately, the 100% efficiency solution is obviously 
rather difficult to excite, since it occurs for only a very 
narrow range of Po values. One approach to dealing with 
this challenge is to relax the requirement of 100% effi- 
ciency (which will never be obtained in practice anyway 
due to losses), and operate at a power Pi < P cr it- In par- 
ticular, Fig.IHJleft) shows that there is a much larger sta- 
ble region for Pi « 0.8P cr it with efficiency around 90%, 
leading one to suspect that this solution may be easier 
to excite than the 100%-efficiency solution at Pi = P cr it- 
This is indeed the case, as is shown in Fig. [TIT right), 
plotting efficiency vs. Po at point C with Pi « 0.8P cr it. 
In this case, there are only two stable solutions, consis- 
tent with Fig. [SJleft), and there are much wider ranges 
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stable) in the phase diagram (Fig. 2}, with the initial con- 
ditions perturbed from the 100%-efficiency stable solution. 
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and 5A 3 o, respectively, with SA^/A?^ = ^30/^. The 
oscillation of the steady-state efficiency with the perturba- 
tion strength is an indication of the complexity of the phase 
space and the shapes of the basins of attraction of each fixed 
point. 



of Po that attain the high-efficiency (~ 90%) solution. 

There are also many other ways to excite the high- 
efficiency solution (or whatever steady-state solution is 
desired). For example, because the cavity is initially de- 
tuned from the input frequency, as described in Sec. HI 



FIG. 10: One way of exciting the system into a controlled 
stable solution: the input power is the sum of an exponen- 
tial turn-on (the blue curve, Pi) and a Gaussian pulse with 
amplitude Po and width ST. The amplitude Po is altered to 
control which stable solution the system ends up in. 



much of the initial pulse power is actually reflected dur- 
ing the transient period, and a more efficient solution 
would vary the pulse frequency in time to match the cav- 
ity frequency as it detunes. One can also, of course, vary 
the initial pulse width or shape, and by optimizing the 
pulse shape one may obtain a more robust solution. 

In particular, one can devise a different (constant- 
frequency) input pulse shape that robustly excites the 
high-efficiency solution, insensitive to small changes 
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appropriate input-pulse amplitude. Right: Same, but for an 
efficiency is ~ 90% from Fig. fright), but is easier to excite. 



in the initial conditions, by examining the power- 
bifurcation diagram in Fig. [Deleft) in more detail. First, 
we observe that input powers > 1.45P cr it have only one 
stable solution, meaning that this stable solution is ex- 
cited regardless of the initial conditions or the man- 
ner in which the input power is turned on. Then, if 
we slowly decrease the input power, the solution must 
"adiabatically" follow this stable solution in the bifur- 
cation diagram until a power w 0.95P cr it is reached, at 
which point that stable solution disappears. In fact, 
by inspection of Fig. [H^left), at that point there are no 
stable solutions, and solution jumps into a limit cycle. 
If the power is further decreased, a high-efficiency sta- 
ble solution reappears and the system must drop into 
this steady state (being the only stable solution at that 
point). This process of gradually decreasing the power 
is depicted in Fig. [T^Heft), where the instantaneous "ef- 
ficiency" is plotted as a function of input power, as the 
input power is slowly decreased. (The efficiency can ex- 
ceed unity,because we are plotting instantaneous output 
vs. input power, and in the limit-cycle self-pulsing solu- 
tion the output power is concentrated into pulses whose 
peak can naturally exceed the average input or output 
power.) Already, this is an attractive way to excite a 
high-efficiency (> 90%) solution, because it is insensitive 
to the precise manner in which we change the power as 
long as it is changed slowly enough — this rate is deter- 
mined by the lifetime of the cavity, and since this lifetime 
is likely to be sub-nanosecond in practice, it is not dif- 
ficult to change the power "slowly" on that timescale. 
However, we can do even better, once we attain this 
high-efficiency steady state, by then increasing the power 
adiabatically. As we increase the power, starting from 
the high-efficiency steady-state solution below the critical 
power, the system first enters limit-cycle solutions when 
the power becomes large enough that the stable solution 



disappears in Fig. [Hf left) . As we increase the power fur- 
ther, however, we observe that these limit cycles always 
converge adiabatically into the 100%-emciency solution 
when P — > P C rit- This process is shown in Fig. [T2lfright). 
What is happening is actually a supercritical Hopf bifur- 
cation at the two points where the upper branch changes 
between stable and unstable: this is a reversible tran- 
sition between a stable solution and a limit cycle (ini- 
tially small oscillations, growing larger and larger away 
from the transition). This is precisely what we observe 
in Fig. (T2J m which the limit cycle amplitudes become 
smaller and smaller as the stable solutions on either side 
of the upper branch are approached, leading to the ob- 
served reversible transitions between the two. The im- 
portant fact is that, in this way, by first decreasing and 
then increasing the power to P cr i t , one always obtains the 
100%-emciency solution regardless of the precise details 
of how the power is varied (as long as it is "slow" on the 
timescale of the cavity lifetime). 



V. CONCLUDING REMARKS 

We have shown that a doubly-resonant cavity not only 
has high-efficiency harmonic conversion solutions for low 
input power, as in our previous work [l|, but also ex- 
hibits a number of other interesting phenomena. We 
showed under what conditions the high-efficiency solu- 
tion is stable, how to compensate for self-phase modu- 
lation, the existence of different regimes of multistable 
solutions and limit cycles controlled by the parameters 
of the system and by the input power, and how to ex- 
cite the desired high-efficiency solution. Although we 
did not observe chaos, it seems possible that this may 
be obtained in future work for other parameter regimes, 
e.g. for pulsed input power as was observed in the 
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case [3]. These dynamical phenomena depend only on 
certain dimensionless quantities a//?, us/u^i, Ts/ri, and 
si+fsi^, although the overall power and time scales de- 
pend upon the dimensionful quantities and so on. 

All of the calculations in this paper were for an ide- 
alized lossless system, as our main intention was to ex- 
amine the fundamental dynamics of these systems rather 
than a specific experimental realization. However, we 
have performed preliminary calculations including both 
linear losses (such as radiation or material absorption) 
and nonlinear two-photon absorption, and we find that 
these losses do not qualitatively change the observed phe- 
nomena. One still obtains multistability, limit cycles, bi- 
furcations, and so on, merely at reduced peak efficiency 
depending on the strength of the losses. In a future 
manuscript, we plan to explore these effects in more de- 
tail in realistic material settings, and propose specific ge- 
ometries to obtain the requisite doubly-resonant cavities. 
In particular, to obtain widely-spaced resonant modes uo 
and 3a; in a nanophotonic (wavelength-scale) context (as 
opposed to macroscopic Fabry-Perot cavities with mir- 
rors), the most promising route seems to be a ring res- 
onator of some sort [42], rather than a photonic crys- 
tal [23[ (since photonic band gaps at widely separated 
frequencies are difficult to obtain in two or three dimen- 
sions). Although such a cavity will naturally support 



more than the two uj\ and cos modes, only two of the 
modes will be properly tuned to achieve the resonance 
condition for strong nonlinear coupling. 

Finally, we should mention that similar phenomena 
should also arise in doubly and triply resonant cavities 
coupled nonlinear ly by sum/difference frequency gener- 
ation (for x^) or four- wave mixing (for x^)- The ad- 
vantage of this is that the coupled frequencies can lie 
closer together, imposing less stringent materials con- 
straints and allowing the cavity to be confined by narrow- 
bandwidth mechanisms such as photonic band gaps [23[ , 
at the cost of a more complicated cavity design. 
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